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Abstract 

This paper describes the procedure for utilizing a damage dependent constitutive 
model to predict progressive damage growth in laminated composites. In this model the 
effects of the internal damage are represented by strain-like second order tensorial damage 
variables and enter the analysis through damage dependent ply level and laminate level 
constitutive equations. The growth of matrix cracks due to fatigue loading is predicted 
by an experimentally based damage evolutionary relationship. This model is incorporated 
into a computer code called FLAMSTR. This code is capable of predicting the constitutive 
response and matrix crack damage accumulation in fatigue loaded laminated composites. 
The structure and usage of FLAMSTR are presented along with sample input and output 
files to assist the code user. 

As an example problem, an analysis of crossply laminates subjected to two stage 
fatigue loading has been conducted herein and the resulting damage accumulation and 
stress redistribution have been examined to determine the effect of variations in fatigue 
load amplitude applied during the first stage of the load history. It is found that the model 
predicts a significant loading history effect on damage evolution. 



Introduction 


Laminated continuous fiber composites are increasingly being utilized in engineering 
applications such as primary load bearing aircraft components. This is in part due to 
the lower weight and higher specific stiffness obtained by using advanced composite ma- 
terials. Unfortunately, laminated composite materials are susceptible to the development 
of microstructural damage when subjected to service loads. This damage includes ma- 
trix cracking, delamination, fiber-matrix interface debonding, and fiber fracture. Each 
microcrack is in itself relatively insignificant since most cracks are arrested at the fibers 
or adjacent plies. However, the resulting redistribution of load to the surrounding re- 
gions creates stress fields favorable to the initiation and propagation of additional damage. 
Catastrophic failure is triggered when the remaining load paths are no longer able to sup- 
port the load. In addition, the strength and stiffness of the material are degraded as a 
result of the load redistribution and the decrease in load paths during the accumulation of 
the subcritical damage. Since the initiation and accumulation of this subcritical damage 
are highly dependent on the stress state within the material, analysis of the structural 
response and service life must account for this stress redistribution. 

Due to the multitude of microcracks often observed in laminated composites, it may 
not be practical to model each flaw explicitly. An alternative approach represents the dis- 
tributed damage by volume averaged quantities known as internal state variables (IS V). 
These quantities may describe the average physical attributes of the distributed damage 
or they might describe the effects of the distributed damage on the material response. In 
the current approach the damaged material volume is modeled as a continuous domain 
with altered properties. Thus, although the microcracks in the representative volume ele- 
ment (RVE) are treated as internal boundaries, the global structural problem is treated as 
simply connected with spatially variable reduced stiffness obtained from a micromechan- 
ics solution for the RVE. As a result, the global problem is made more computationally 
tractable. Even though the ply level stresses obtained from this approach are locally av- 
eraged quantities, the model predictions are in qualitative agreement with experimental 
results. This approach, call damage mechanics, is suitable for damage that is small in size 
relative to the scale of the structure being analyzed and is spatially homogeneous within 
the RVE. Although it has been applied to the study of a wide range of phenomena from 
microcracking to chemical and radiation damage in engineering materials, only recently 
has it been applied to matrix cracking in laminated composites [1-9]. 

This report describes how to perform an analysis using this approach. The model is 
formulated for matrix cracking and delamination damage in laminated polymeric compos- 
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ites with a brittle matrix. At the present time, a damage evolutionary relationship for 
fatigue induced matrix cracking has been implemented into the lamination computer code, 
FLAMSTR. This code is used to demonstrate some of the effects that the loading history 
has on the accumulation of matrix cracking. The evolutionary relationships for delamina- 
tion damage are currently under development and will be incorporated into the code at a 
future time. A brief description of this model will be presented. Detailed development of 
this model can be found in the published literature [5,6,10-14], 

Model Description 

In the proposed model, the effects of the matrix cracks are introduced into the ply 
level constitutive equations as follows [15]: 

Wl.j} = [QKfiLy - <*£(,} t 1 ) 


where cr i are the locally averaged components of stress, [Q] is the ply level transformed 
stiffness matrix, and e/ MJ are the locally averaged components of strain. a£* are the 
components of the strain-like internal state variable for matrix cracking and are defined 
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where V'/. is the volume of an arbitrarily chosen representative volume of ply thickness which 
is sufficiently large that atj* do not depend on Vl, are the crack opening displacements, 
and n j are the components of the vector normal to the crack face. It is assumed in 
the current model that the internal state variable representing the mode I matrix 

crack opening, is the only nonzero component. For a uniaxially loaded medium containing 
alternating 0° and 90° plies, a^, has been found from a micromechanics solution to be 
related to the far field normal force and crack spacing as follows [15]: 
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where 

x no j 

p is the force per unit length applied normal to the fibers and 2 i and 2a are the layer 
thickness and crack spacing, respectively. C 1 2222 the modulus in the direction transverse 

to the fibers and C 1212 * s the inplane shear modulus. Both moduli are the undamaged 
properties. 


3 



Equation (3) requires that the matrix crack spacing be known in each ply of the 
laminate. Since it is usually necessary to predict the damage accumulation and response 
for a given loading history, damage evolutionary relationships must be utilized to determine 
the values of the internal state variables. The authors have used the following relationship 
for the rate of change of the internal state variable in each ply during fatigue loading 
conditions and when the available strain energy release rate is greater than some critical 
value (? c , [16], 


da 
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22 
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(5) 


da M 

where describes the change in the internal state variable for a given change in the 

crack surface areas, k and h are material parameters, and .V is the number of load cycles. 
G is the damage dependent strain energy release rate for the ply of interest and is calculated 
from the following equation, 


G = ViCijki{sL lJ 
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where Vi is the local volume. Interactions with the adjacent plies will result in ply strains, 
sl which are affected by the strains in adjacent plies. Thus, the energy release rate, G, 
in each ply will be implicitly reflected in the calculation of the ply level response, so that 
equation (5) is not restricted to a particular laminate stacking sequence. Utilizing equation 
(6) in equation (5) and integrating the result in each ply over time thus gives the current 
damage state in each ply for any fatigue load history. 

The ply level strains are defined as follows: 
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where e° L and kl are the midplane strains and curvatures, respectively. The aforementioned 
ply strains are then substituted into equation (1) to produce the ply level stresses. 
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Damage dependent lamination equations are obtained by integrating these ply stresses 
through the thickness of the laminate jll|. Next, the stiffness matrix in the laminate 
equation is inverted to produce: 


where [.-I], 1 B], and \D\ are, respectively, the undamaged laminate extensional, coupling, 
and bending stiffness matrices. They are defined by [17j, 

Ml “ [Q\k( z k - -*-i) ( 14 ) 
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where [Qj*. is the elastic modulus matrix for the A-" 1 ply in laminate coordinates. N are 
the components of the resultant force per unit length and M are the components of the 
resultant moments per unit length; {/ u } and {gr' u } represent the contribution to the 
resultant forces and moments from matrix cracking and are calculated from, 

( l7) 
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where {<v' u )n. contains the matrix cracking internal state variables for the k th ply. Thus 
given the forces, N, and moments, M, as well as the damage variables in each ply, equation 
(13) can be utilized to calculate the midsurface strains, and curvature, 

Program Description 

Program Structure 

The damage dependent lamination model has been coded into the FORTRAN program 
FLAMSTR following the algorithm shown in Figure 1. This program enables the analysis 
of the stress-strain response and accumulation of matrix cracking at a material point in 
a laminate subjected to fatigue loading. The program begins by reading in the laminate 
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Figure I. Structure of FLAMSTR used in the analysis of damage accumulation. 
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and material properties from a data file. The extensional, coupling, and bending stiffness 
matrices, [ t|, [B], and [£)], respectively are calculated and then inverted for later use. 
Next, the initial damage states for each ply and the growth law parameters are read into 
the program. They can be determined if the initial crack spacing is known from equation 
(3). The procedure for the determination of the damage evolution parameters can be 
found in Appendix A. Finally, the fatigue loading history is read in. The constitutive 
portion of the program initiates with the calculation of the “damage induced forces and 
moments using equations (17) and (18), respectively. These quantities are combined with 
the applied forces and moments to calculate the midplane strains and curvatures through 
equation (13). The strains in each ply are then determined with equations (7) through 
(12). Equation (1) is then used to calculate the damage dependent stresses in each ply 
and the strain energy release rate is determined for each ply from equation (6) using the 
ply stresses. The evolution of a™ in each ply is then obtained from equation (5) and the 
damage state is updated. The above procedure is repeated for the desired number of load 
history increments. 

The computational algorithm assumes that the rate of damage evolution is small 
enough that the strain energy release rate can be considered to be constant during each 
load cycle. Often this condition leads to exponential overflow errors during the execution 
of the computer code. The high sensitivity of the mode I matrix cracking ISV to the strain 
energy release rate is due to the power law form of the damage evolutionary relationship. 
To model this change in the strain energy release rate with the current algorithm, the load 
cycles experiencing large changes in damage evolution are divided into subincrements for 
calculation. This enables the strain energy release rate to reflect the damage accumulated 
during the load cycle. This approach has produced satisfactory results. Subincrementation 
has also been found to be necessary during the initial cycles of the fatigue load history 
when the laminate goes from an undamaged to a damaged state and whenever the fatigue 
load amplitude is increased. 

Program Inputs 

The execution of FLAMSTR requires the creation of two input data files. One of these 
files is labeled datfil.d and contains two entries. The first entry is the name of the other 
input file from which all the information required in the calculations are retrieved. On the 
next line is the name of the data file in which the output will be stored. This output file 
is created by the program during execution. The second input file contains information 
about the stacking sequence, material properties, initial damage state, and the loading 
history. The variables, listed in the order in which they should appear in the data file, 
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along with their description are presented in the following list: 

Number of plies in laminate. 

Components of the transformed ply stiffness matrix. 
Components of the transformed ply stiffness matrix. 

Damage condition: 0 — no damage, 1 = matrix cracks. 
Applied forces. 

Applied moments. 

/(/), theta(i) Thickness and orientation (deg) of ply i. For i = 1, nplies. 

alphnm(i ,2), alpharn(i, 8) Initial values of the Mode I and Mode II matrix cracking ISV 

for ply i. For i - 1, nplies. 

Slope of the relationship between the far field normal stress 
and and growth law parameters. 

Initial cycle number; final cycle number; number of increments 
taken to go from nci to nr/; I ■ nine ( nef — nci ). 

Output results to datafile every iprnum increments. Number 
of subincrements employed during a change in the maximum 
load. 

Cycle number at which the maximum applied load is changed; 
maximum applied load factor. ( nci njump < nef, set xfac 
— 1 if applied load remains constant ) 

Note that each line is read in a format free manner. A sample set of input files can be 
found in Appendix B. 

Program Outputs 

The output data file contains a listing of the input variables plus the program gen- 
erated results. Due to the number of load history increments involved in the analysis 
of fatigue, only the results at preselected increments are stored in the output file. The 
frequency at which the data is recorded is set by the variable iprnum in the input files. In- 
formation stored consists of the load cycle number, the values of the matrix cracking ISV s 
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in each ply, the midplane strains and curvatures, and the damage dependent strains and 
stresses at the outer fiber of each ply. The output file corresponding to the aforementioned 
sample input files can also be found in Appendix B. 

Sample Calculations 

The model has been employed to simulate the damage dependent stress redistribution 
and the accumulation of matrix cracks in a {0/902 /0J a laminate subjected to the uniaxial 
loading histories shown in Figure 2. Each load case consists of two load segments with 
the latter at a maximum load of 600 lb/in. A fatigue load ratio of 0.1 has been assumed 
for both load segments. Material properties for AS4/3501-6 Graphite/Epoxy, as listed in 
Table l, have been used in the calculations. The sample input and output FLAMSTR files 
found in Appendix B correspond to load case I. 

Figure 3 illustrates the dependence of the accumulation of matrix cracks on the max- 
imum load amplitude. The damage caused by loading initially at a lower load amplitude, 
cases I and 11, causes only a minor effect on the final damage state as compared to load case 
111 in which a maximum load of 600 lb/in is maintained during the entire load history. 
However, when the initial fatigue load level is 800 lb/in (case IV ) much of the damage 
occurs during the initial fatigue cycles. Furthermore, even though only a minute amount 
of damage accumulates during the second load segment, the amount of damage at the 
completion of this fatigue load history is almost twice as much as the other three cases. 
The decelerated growth of matrix cracks during the second stage occurs when the spacing 
between the matrix cracks no longer enables the transfer through shear of sufficient load 
back to the ply to create additional matrix cracking. On the other hand, the negligible 
accumulation of damage during the first stage of load in case I, when the maximum load is 
200 Ib/in, indicates that the applied load is not sufficient to produce an appreciable amount 
of damage. Thus, even though a critical value for damage growth has not been specified 
in the calculations, the damage evolution relationship shown by equation (5) behaves as if 
a threshold exists for damage growth. The accuracy of this predicted behavior, however, 
requires further investigation. 

The average axial stresses in the 90° plies of the crossply laminate are shown in 
Figure 4. The plies experience a sharp decrease in axial stress upon the initiation of 
matrix cracking. This is followed by a period of gradual decrease as further damage 
accumulates. The axial ply stress then abruptly changes when the maximum load is 
changed to 600 lb/in at the initiation of the second load segment. Since each case has 
the same maximum fatigue load the corresponding axial stress serves as an indication of 
the relative load carrying capability remaining in the 90° plies. The laminate retains a 
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Figure 2. Fatigue load amplitude applied to crossply laminates. 




Table 1. Material Properties for AS4/3501-6 
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Figure 4. Axial stress in the 90” plies. 
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similar amount of load carrying capability for the first three load cases. However, the load 
carrying capability for load case IV is substantially reduced compared to the other cases. 
This is in accordance with the amount of damage sustained for this load case. 

The results presented illustrate the damage accumulation and stress redistribution oc- 
curring within uniaxially fatigue loaded crossply laminates. As the load carrying capability 
diminishes in the 90° plies, the surrounding 0° plies must assume a greater portion of the 
applied load. The load transfer continues until the axial stress in the 0 ° plies exceeds the 
strength of a fiber. Since the 0° plies are the primary load bearing member for crossply 
layups, laminate failure is probably triggered by this event. Although this example is a 
relatively simplistic representation of the progressive failure process found in laminated 
composites, it does demonstrate how this model can be employed to analyze this process. 
The ability to model the damage dependent stress redistribution and damage accumulation 
can be advantageous in the design and maintenance of laminated composites. By simulat- 
ing the progressive failure of a laminate, potential damage modes and their locations can 
be identified. The laminate can thus be redesigned to suppress or minimize the effects of 
such damage. 

The introduction of multiaxial loading and/or angle plies will promote the develop- 
ment of mixed mode matrix cracking and delamination during the application of fatigue 
loads. In addition, fiber breakage in one ply may not cause the unstable fracture of the 
laminate. While this model possesses the capability to account for mixed mode matrix 
cracking and delamination, relationships describing their evolution are necessary to analyze 
their contribution in the progressive failure process. The present version of FLAMSTR 
calculates only the mode I matrix crack contribution in each ply. This factor must be con- 
sidered in the interpretation of the predicted results until this code is updated to account 
for these conditions. To analyze laminated composite structures with spatially varying 
stress fields this constitutive code has been incorporated into a finite element structural 
analysis algorithm [18]. This finite element analysis code has been successfully used to 
predict spatial variations of damage evolution in structural components such as a plate 
with a circular cutout [19]. 


Summary 

A computer code called FLAMSTR has been developed for the analysis of the pro- 
gressive failure process in fatigue loaded laminated composites. This code utilizes volume 
averaged internal state variables to describe the kinematics of matrix cracks. The effects 
of the damage enter the analysis via the damage dependent laminate equations. Ply level 
stresses obtained from the results of these equations are then used to predict the amount 
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of mode I matrix crack growth. This procedure is repeated for each load cycle to simulate 
the damage accumulation process. 

This report serves as a user’s guide for the constitutive code. Sample input and output 
datafiles are enclosed to assist in the usage of this computer code. A simulation of crossply 
laminates subjected to two stage uniaxial fatigue load histories has been conducted to 
illustrate the ability of the model to predict the path dependent damage accumulation and 
stress redistribution in the laminate. Since the information generated from this analysis 
can be used to construct the sequence of events leading to the failure of the laminate, it 
has potential uses in the design of composite structures. 
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Appendix A. Determination of Growth Law Parameters 

j M 

The damage evolution parameters, ^c 2 , k, and ft, are unique to each material system 

da M 

and must be determined prior to the analysis. The term -j^ 2 reflects the changes to 
the internal state variable with respect to changes in the crack surfaces. If it is assumed 
that the matrix crack surfaces are flat and aligned perpendicular to the plane formed by 
the ply, thereby permitting the description of the total crack surface area in terms of the 

, M 

crack spacing, then can be calculated analytically form equation (3). It was found, 
for typical brittle graphite/epoxy material properties, that 022 exhibits an almost linear 
relationship to the crack spacing at each far field load. Thus this analysis assumes the 

slope, '-jo 2 ) to be constant for each far field load. Furthermore, equation (3) indicates that 

I XI . , dct M 

-j|2 varies linearly with the far field load. It is the slope of the relationship between -jj?- 
and the far field load that is required by FLAMSTR. This variable is stored in the input 
file under the named dpara. 

The constants k and ft in equation (5) must be determined from experimental data. 
Since k and ft are assumed to be parameters unique to each material system, the values 
determined from one laminate stacking sequence and one loading condition should be valid 
for all other cases. This has been shown for uniaxially fatigue loaded crossply laminates 
with varying number of consecutive transverse plies and maximum fatigue stresses [20]. 

To evaluate k and ft, a curve of the quantity versus the strain energy release 

rate, G, must be generated for a particular stacking sequence and maximum fatigue stress. 
Experimental data from uniaxially fatigue loaded crossply laminate are used because it can 
be assumed that the matrix crack opening mode to be essentially mode I. Thus the ISV, 
<* 22 , > s sufficient in describing the damage state and the contributions from the other two 
opening mode do not have to be considered in the calculations. In addition, the transverse 
plies should be grouped together in the laminate so that the damage evolution at a single 
layer needs to be considered. Since most damage accumulation data are reported as crack 
spacing or density at a particular point in the loading history, equation (3) is utilized to 
convert this data to the form of the ISV, <*22 • The resulting <*22 versus load cycle, N, curve 
serves as the starting point in this procedure. The following steps describe this process: 

(1) <I jn l at a point in the loading history is determined by taking the slope of the 

vs. N curve. This task can be facilitated if the 0*2 vs. N curve is fitted numerically 
and the first derivative taken. 

j J VI 

(2) Equation (3) is employed to calculate -jj?- at each data point. These will be used in 
the determination of the strain energy release rate and growth law parameters. 

(3) The strain energy release rate, G , is calculated using equation (6) where the height of 
the local volume, Vj, is assumed to be equal to the thickness of the damaged layer. 
Since and the maximum applied fatigue stress are known, si tJ in equation (6) can 
be determined via equations (7) through (12). This is performed for each data point. 

(4) Now -j& is divided by the corresponding so that the damage evolution law, 
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equation (9), becomes 


da 22 dS 

Tn 

With the left hand side of equation (Al) and the corresponding values of the strain 
energy release rate, G, calculated at various points in the fatigue loading history, k and 
n can be determined by taking the natural log of equation (Al) and then employing 
a linear regression procedure. 



For the present model, k and n of AS4/3501-6 graphite/epoxy are determined from 
damage accumulation data published by Chou, et al. [21]. The data employed is for 
[0 2 /90 2 ]« AS4/3501-6 laminates fatigue loaded at the maximum fatigue stress of 43 ksi 
and a stress ratio of 0.1 . The parameters are determined to be 

k = 4.42, h = 6.39 . (A2) 

This values are then used in the calculations involving the AS4/3501-6 graphite/epoxy 
system. 
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Appendix B. Sample Input and Output Files 


Input datafile: datfil.d* 

finp.dat iname 

ft200.d outname 


Input datafile: finp.dat* 

6 


20.29eG 

0.433eG l.GleG 

O.433e0 

l.GleG 0.695eG 

1 


200.0 0.0 0.0 

0.0 0.0 

0.0 

0.00505 

0.0 

0.01010 

90.0 

0.00505 

0.0 

0.00505 

0.0 

0.01010 

90.0 

0.00505 

0.0 

0.0 0.0 


0.0 0.0 


0.0 0.0 


0.0 0.0 


0.0 0.0 


0.0 0.0 


6.9017e-G 

4.41613 6.388592 

0.0 10000.0 10000.0 

500 200.0 

5000.0 3.000 


n plies 

Q I 1 j Q 1 2 i Q 22 
Q 13 , Qsi > Q 66 

if lag 

N x , N y , N iy 
iV/j. j M y , M ry 
f(l), tluta(l) 
t{ 2), //>cfa(2) 

/( 3 ), thttaf 3 ) 

/(4), f/iefa(4) 

/(5), theia.(5) 
f(6), theta{ 6) 

« /p/» a m ( 1 , 2 ) , alph a m ( 1 , 8 ) 
alpham( 2, 2), alpham( 2, 8) 
alpham( 3, 2), alpham(3, 8) 
alpham{A, 2), alphamf 4,8) 
alpham(5, 2), alpham(5,8) 
alpham(6, 2), alpham(6, 8) 
dpara, k, h 
nci, ncf , mnc 
iprnum, nsubic 
n jump, xfac 


* Italicized items are descriptions of the input variables and are not part of the data 
file. 
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Output datafile: ft200.d 



OUTPUT FOR 


6 PLIES 


THE C MATRIX IS 


0. 2029000E+08 0 . 4330000E +06 

0. 4330000E+06 O. 1610000E+07 


0. 1 6 10000E+07 
0. 6950000E+06 


I FLAG * 1 

THE APPLIED FORCE Ni ARE 

0 . 2 OOOOOO E + 0 3 O . OOOOOOOE+OO 0 . OOOOOOOE+OO 0 . OOOOOOOE+OO 
0 . OOOOOOOE+OO 0 . OOOOOOOE +00 


THE APPLIED MOMENT Mi ARE 

0 . OOOOOOOE+OO 0 . OOOOOOOE+OO 0 . OOOOOOOE + OO 0 . OOOOOOOE +00 
0 . OOOOOOOE +00 0 . OOOOOOOE+OO 


PLY NO. 


T 


THETA 


1 

2 

3 

4 

5 

6 


0.5050000E-02 
0. 1010000E-01 
0. 5050000E-02 
0 . 5050000E -02 
0. 1010000E-01 
0. 5050000 E- 02 


0. OOOOOOOE +00 
0 . 9000000 E + 02 
0. OOOOOOOE+OO 
0. OOOOOOOE +00 
0 . 9000000E+02 
O. OOOOOOOE +00 


DPARA s 0 . 690 1 700E -05 XK 1 * 0 . 44 16 1 30E+0 1 XN1* 0 . 6388592 E+0 1 


INITIAL CYCLE: 0 . OOOOOOOE +00 FINAL CYCLE: 0 . 1000000E+05 

CYCLE INCREMENT: 0 . 1 OOOOOOE+O 1 

PRINT OUTPUT INCREMENT: 500 

SUBINCREMENT CYCLES (RAMP UP): 0 . 2000000E+03 


LOAD JUMP AT CYCLE: 0 . 5000000E +04 LOAD FACTOR: 


0. 3000000E +0 1 


CYCLE NUMBER 0 . 4990000E+03 


PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0. 182843 1 D - 18 
0.7381833D-08 
0. 182843 ID- 1 8 
0. 1 82843 1 D — 1 8 
0 . 738 1833D-08 
0. 1 8284 3 1 D- 1 8 


ALPHAM8 

0 . OOOOOOOD + 00 
-O. 32267020 - 1 5 
0 . OOOOOOOD +00 
0 . OOOOOOOD +00 
-0 . 3 2 26 702 D- 1 5 
0 . OOOOOOOD +00 


EPSO ( 1 ) - EPSO( 6 ) 


0.4539266D-03 
0 OOOOOOOD +00 

KAPPA ( 1 ) - KAPPA ( 6 ) 

-0 . 295 1 98 1 D- 1 7 
0 . OOOOOOOD +00 


*0. 1 67874 1 D“04 
0. OOOOOOOD +00 


0 . 1 0700 1 1 D - 1 8 
0 . OOOOOOOD +0O 


-0 . 587830 ID -04 
-0.6708493D- 1 1 


0. 1293403D- 18 
0 . 480 1 262D-25 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E(X) E(Y) E(Z) E(XY) 


1 

2 

3 

4 


0.45393D-03 

0.45393D-03 

0-45393D-03 

0.45393D-03 


-O. 16787D-04 
-0. 16787D-04 
-0. 16787D-04 
-0. 16787D-04 


-0.58783D-04 
-0.58783D-04 
-0. 58783D-04 
-0.58783D-04 


-0.67085D-1 1 
-0.67085D-1 1 
-O . 67085D- 1 1 
-0 . 67085D- 1 1 
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5 0 45393D-03 -0.167870-04 -0.587830-04 -0670850-11 

6 0. 453930-03 -0.16787D-04 -0.58783D-04 -0.67085D-11 

THE STRESSES AT THE OUTER FIBER ARE 


PLY 

NO. 

S( 1 ) 

S ( 2 ) 

$(3) 

S(( 


1 

9177.45 

169.52 

101.91 

0.00 


2 

-169.52 

723.54 

-101 .91 

0.00 


3 

9177.45 

169.52 

101.91 

0.00 


4 

9177.45 

169.52 

101.91 

0.00 


5 

-169.52 

723 . 54 

-101.91 

0.00 


6 

9177.45 

169 . 52 

101.91 

0.00 

CYCLE NUMBER 0. 

9990000E+03 



PLY 

NO. 

ALPHAM2 

ALPHAM8 



1 

0.36104510- IS 

0 . OOOOOOOD+OO 



2 

0. 14574560-07 

-0. 6370742D- 15 



3 

0. 36104510-18 

0 . OOOOOOOD+OO 



4 

0.36104510-18 

0 . OOOOOOOD+OO 



5 

0. 14574560-07 

-0 . 6370742D- 15 



6 

0.36104510-18 

0. OOOOOOOD+OO 



EPSO ( 1 ) - EPSO( 6 ) 

0.45392710-03 
0 . 00000000+00 

KAPPA ( 1 ) - KAPP A ( 6 ) 

-0 . 2952008D- 1 7 
0 - OOOOOOOD+OO 


-0. 1678728D-04 
0 . OOOOOOOD+OO 


O. 10699160- 18 
0 . 00000000+00 


-0. 58783090-04 
-0.67084180-1 1 


0. 1 293454D- 1 8 
0 . 480080 1 D - 2 5 


THE STRA I NS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E(Y) E(Z) E(XY) 


1 

2 

3 

4 

5 

6 


0.453930-03 
0.453930-03 
0 . 453930-03 
0. 45393D-03 
0 45393D-03 
0. 45393D-03 


-0. 16787D-04 
-0. 16787D-04 
-0. 167870-04 
-O 167870-04 
-O. 167870-04 
-0. 167870-04 


-0.587830-04 

-0.587830-04 

-0.587830-04 

-0.587830-04 

-0.587830-04 

-0.587830-04 


-0.670840- 1 1 
-0.670840-1 1 
-0.67084D-1 1 
-0.670840-1 1 
-0.67084D-1 1 
-0.67084D-1 1 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO. S< 1 ) S ( 2 ) S ( 3 ) S(6) 


1 

9177 . 46 

169.52 

101.91 

0.00 

2 

-169.52 

723.53 

-101.91 

0.00 

3 

9177.46 

169.52 

101.91 

0.00 

4 

9177.46 

169.52 

101.91 

0.00 

5 

-169.52 

723.53 

-101.91 

0.00 

6 

9177.46 

169.52 

101.91 

0.00 


CYCLE NUMBER 0 . 1 499000E+04 


PLY NO. 

ALPHAM2 

ALPHAM8 

1 

0.53925330-18 

0 . 00000000+00 

2 

0. 2176582D-07 

-0.9514 1 42D- 1 5 

3 

0.53925330-18 

0 . OOOOOOOD+OO 

4 

0. 53925330-18 

0 . OOOOOOOD+OO 

5 

0. 21765820-07 

-0.95141 4 2D- 1 5 

6 

0.53925330-18 

0 , OOOOOOOD+OO 

EPSO(I) - 

EPSO ( 6 ) 



0.45392770-03 
0 . 00000000+00 


-0. 16787160-04 
0 . 00000000+00 


-0.5878318D-04 
-0.6708343D-1 1 


KAPPA ( 1 ) - KAPPA ( 6 ) 

-0. 2952028D- 17 
0. OOOOOOOD+OO 


0. 1069888D- 18 
0 . OOOOOOOD+OO 


0. 12934880-18 
0.48006580-25 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E ( Y ) E(Z) E(XY) 


1 

2 

3 

4 


0.453930-03 
0.45393D-03 
0. 453930-03 
0.453930-03 


-0. 16787D-04 
-0. 167870-04 
-0. 167870-04 
-0. 16787D-04 


-0. 587830-04 
-0.587830-04 
-0 - 587830-04 
-0.587830-04 


-0.67083D- 1 1 
-0.67083D- 1 1 
-0.6708 3D- 11 
-0670830-11 
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5 O.45393D-03 -0.167870-04 -0.58783D-04 

6 O . 45393D-03 -O . 1 6787D -04 -0.58783D-04 


THE STRESSES AT THE OUTER FIBER ARE 


PLY NO 

S( 1 ) 

S ( 2 ) 

SC 

1 

9177.49 

169.52 

101.91 

2 

-169.52 

723.50 

-101 .91 

3 

9177 .49 

169.52 

101.91 

4 

9177.49 

169 . 52 

101.91 

5 

- 169 . 52 

723.50 

-101 .91 

6 

9177.49 

169.52 

101.91 

CYCLE 

NUMBER 0. 

2999000E+04 



PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0. 1073915D- 17 
0. 43330840-07 
0. 10739150-17 
0. 10739150-17 
0. 4333084D-07 
0. 10739150-17 


ALPHAM8 

O.OOOOOOOD+OO 
-0. 18940510- 14 
O . 00000000+00 
O.OOOOOOOD+OO 
-0. 1 89405 1 D - 1 4 
0 . 00000000+00 


EPSO ( 1 ) - EPSO ( 6 ) 


0.45392920-03 -0 . 1 678679D -04 -0 

0 . OOOOOOOD +00 0 . 00000000 +00 - O 


KAPPA ( 1 ) - KAPPA ( 6 ) 

-0.2952133D-17 
0 . 00000000+00 


0. 1069506D- 18 0 

O.OOOOOOOD+OO 0 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E(Y) E(Z) 


1 

2 

3 

4 

5 

6 


0.45393D-03 
0 • 45393D -03 
0.45393D-03 
0.45393D-03 
0.45393D-03 
0.453930-03 


-O. 167870-04 
-0. 167870-04 
-0. 16787D-04 
-0. 167870-04 
-0. 167870-04 
-0. 167870-04 


-0- 58783D-04 
-0. 587830-04 
-0.587830-04 
-0.58783D-04 
-0.587830-04 
-0. 587830-04 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO. S ( 1 ) 5(2) S ( 3 ) 


1 9177.50 

2 -169.52 

3 9177.50 

4 9177.50 

5 -169.52 

6 9177.50 


169.52 
723.49 

169 . 52 

169 . 52 
723.49 

169.52 


101.91 
-101 .91 

101.91 

101.91 
-101 .91 

101.91 


CYCLE NUMBER 0 . 3499000E+04 


PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0. 1252148D-17 
0.50516250-07 
0- 12521480-17 
0. 1252148D-17 
0. 5051625D-07 
0. 1252148D- 17 


ALPHAM8 

O.OOOOOOOD+OO 
-0. 22081360-14 
O.OOOOOOOD+OO 
0 . OOOOOOOD +00 
-0. 2208136D- 14 
O.OOOOOOOD+OO 


EPSO ( 1 ) - EPSO ( 6 ) 

0. 4539298D-03 
0 . 00000000+00 

KAPPA ( 1 ) - KAPPA ( 6 ) 

-0.29522180-17 

O.OOOOOOOD+OO 


-0. 1678667D-04 
0 . 0000000 D + 00 


0. 10692600-18 
0 . OOOOOOOD +00 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E ( Y ) E(Z) 


1 

2 

3 

4 


0.45393D-03 
0453930-03 
0. 45393D-03 
0. 45393D-03 


-O. 167870-04 
-O. 167870-04 
-O. 167870-04 
-O . 167870-04 


-0.587840-04 

-0.587840-04 

-0.58784D-04 

-0.587840-04 


-0.67082D-1 1 
-0.670820-1 1 


S ( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


5878344D-04 
67081 180-1 1 


1 293694D- 18 
47988030-25 


E(XY) 

-0.670810- 1 1 
-0.67081D-1 1 
-0.67081D-1 1 
-0 . 6708 ID- 1 1 
-0.670810-1 1 
-0 . 6708 ID- 1 1 


S ( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


. 5878353D-04 
. 6708043D- 1 1 


. 1 293855D- 1 8 
. 4797596D-25 


E ( X Y ) 

-0 . 67080D- 1 1 
-0 . 67080D- 1 1 
-0.67080D-1 1 
-0.67080D-1 1 
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5 

0. 453930-03 

-0. 16787D-04 

-0.58784D-04 

-0.670800-1 1 

6 0 . 45393D-03 

STRESSES AT THE 

-0. 16787D-04 -0.58784D-04 

OUTER FIBER ARE 

-0.67080D-1 1 

NO. 

S( 1 ) 

S ( 2 ) 

S ( 3 ) 

5(6) 

1 

9177.51 

169 . 53 

101.91 

0.00 

2 

-169.53 

723.48 

-101.91 

0.00 

3 

9177.51 

169.53 

101,91 

0.00 

4 

9177.51 

169 . 53 

101.91 

0.00 

5 

-169.53 

723 48 

-101.91 

0.00 

6 

9177.51 

169 . 53 

10191 

0.00 


CYCLE NUMBER 


0. 3999000E+04 


PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0. 1430387D- 17 
0.577002 ID-07 
0. 1430387D- 17 
0. 1430387D-17 
0.577002 ID-07 
0. 1 430387D- 1 7 


AL PHAM 8 

0 . OOOOOOOD + 00 
-0. 2522 156D- 14 
0 . OOOOOOOD +00 
0 . OOOOOOOD +00 
-0. 2522 156D- 14 
0 . OOOOOOOD *00 


EPSO ( 1 ) - EPSO ( 6 ) 

0. 4539303D -03 
O . OOOOOOOD +00 


-0 1678655D-04 
0. OOOOOOOD +00 


-0758783520-04 
-0.6707968D-1 1 


KAPPA( 1 ) - KAPP A ( 6 ) 

-0. 2952227D- 17 
0 . OOOOOOOD +00 


O. 1069154D- 18 
0 . OOOOOOOD +00 


0. 1293880D- 18 
0 . 4797095D- 25 


THE STRAINS AT THE OUTER FIBER ARE 
PLY NO. E ( X ) E ( Y ) 


O . 45393D-03 
0.45393D-03 
0.45393D-03 
0.45393D-03 
0 . 45393D -03 
0.453930-03 


-0 

-0 

-0 

-0 

-0 

-o 


16787D-04 

16787D-04 

16787D-04 

16787D-04 

16787D-04 

16787D-04 


E( 2 ) 


-0. 58784D-04 
-0.58784D-04 
-0.58784D-04 
-0.58784D-04 
-0 . 58784D-04 
-0.58784D-04 


E(XY) 

-0.67080D- 1 1 
-0.67080D-1 1 
-0.67080D-1 1 
-0 . €70800- 1 1 
-0.67080D-1 1 
-0.67080D-1 1 


STRESSES AT THE 
NO . S ( 1 ) 

OUTER FIBER 
S ( 2 ) 

ARE 

S ( 3 ) 

5 C € 

1 9177.52 

. 169.53 

101.91 

0.00 

2 -169.53 

723 . 47 

-101.91 

0.00 

3 9177.52 

169.53 

101.91 

0.00 

4 9177.52 

169.53 

101.91 

0.00 

5 -169.53 

723 . 47 

-101.91 

0.00 

6 9177.52 

169.53 

101.91 

0.00 


CYCLE NUMBER 


0. 4499000E+04 


PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0. 1608632D- 17 
0 . 648827 ID-07 
0. 1 608632D- 1 7 
0. 1608632D- 17 
0 . 648827 ID-07 
0. 1608632D- 17 


ALPHAM8 

0. OOOOOOOD +00 
-0.28361 1 3D- 14 
0 . OOOOOOOD +00 
0. OOOOOOOD +00 
-0.28361 13D-14 
0 . OOOOOOOD+OO 


EPSO( 1 ) - EPSO( 6 ) 

0. 4539308D-03 
0 . OOOOOOOD+OO 


-0. 1678642D-04 
0 . 00000000+00 


-G.587B370D-04 
-0.6707893D-1 1 


KAPPA ( 1 ) - KAPPA ( 6 ) 

-0 . 2952305D- 1 7 
0. OOOOOOOD+OO 


O. 1069002D- 18 
0 . OOOOOOOD+OO 


0. 1 294020D- 1 8 
0.4796336D-25 


THE 

PLY 


STRAINS AT 
NO . 


THE OUTER 

e ( x ) 


FIBER ARE 
E( Y) 


1 0.45393D-03 

2 0.45393D-03 

3 0.45393D-03 

4 0.45393D-03 


-O. 167860-04 
-O. 167860-04 
-O. 16786D-04 
-O. 16786D-04 


E ( Z > 

0587840-04 
0 . 58784D -04 
0.58784D-04 
0.58784D-04 


E ( X Y ) 

-0.67079D-1 1 
-0.67079D- 1 1 
-0.67079D-1 1 
-0.67079D-1 1 
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5 

0. 45393D-03 

-0. 16786D-04 

-0. 58784D-04 


6 

0. 45393D-03 

-0. 1G786D-04 

-0. 587840-04 

THE 

STRESSES AT THE 

OUTER FIBER ARE 


PLY 

NO 

S( 1 ) 

S ( 2 ) 

5(3) 


1 

9177.53 

169.53 

101.91 


2 

- 169 . 53 

723.46 

101.91 


3 

9177.53 

169 53 

101.91 


4 

9177.53 

169.53 

101.91 


5 

-169.53 

723 46 

101.91 


6 

9177.53 

169 . 53 

101.91 

CYCLE 

NUMBER 

0 . 4999000E+04 


PLY 

NO 

ALPHAM2 ALPHAM8 

1 0. 1786883D- 17 0 . OOOOOOOD+OO 


2 0. 7206375D-07 

3 0. 1786883D- 17 

4 0. 1786883D-17 

5 0.72063750-07 

6 0. 1786883D-17 


-0. 31500070-14 
0 . OOOOOOOD+OO 
0 . 00000000 + 00 
-0.3 1500070-14 
0. OOOOOOOD+OO 


EPSO( 1 ) - EPS0(6 ) 

0. 4539313D-03 
0 . OOOOOOOD+OO 


-0. 16786300-04 
0 . OOOOOOOD+OO 


-O 

-0 


KAPPA ( 1 ) - KAPPA ( 6 ) 

- 0 . 29524280 - 17 0 . 10687080-18 0 

0 . OOOOOOOD+OO O . OOOOOOOD+OO 0 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO, E ( X ) E ( Y ) E(Z) 


1 

2 

3 

4 

5 

6 


0 453930-03 
O. 45393D-03 
0.45393D-03 
0.453930-03 
O. 45393D-03 
0.45393D-03 


-0. 16786D-04 
-O. 16786D-04 
-0. 16786D-04 
-0. 16786D-04 
-0. 167860-04 
-0. 16786D-04 


-0. 58784D-04 
-0.58784D-04 
-0.58784D-04 
-0.58784D-04 
-0.58784D-04 
-0. 58784D-04 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO. S ( 1 ) S(2) 5(3) 


1 9177.55 

2 -169.53 

3 9177.55 

4 9177.55 

5 -169 53 

6 9177.55 


169.53 
723 . 44 
169.53 
169 . 53 
723 . 44 
169 . 53 


101 .91 
-101.91 
101.91 
101.91 
-101.91 
101.91 


CYCLE NUMBER 0 . 5499000E+04 


PLY NO. 
1 
2 

3 

4 

5 

6 


ALPHAM2 
0. 13138560-1 1 
0.50862560-03 
0. 13138560-1 1 
0. 1313856D- 1 1 
0. 5086256D-03 
0. 13138560- 1 1 


ALPHAM8 

0 . 00000000+00 
-0. 2223273D- 10 
0 . OOOOOOOD+OO 
0 . OOOOOOOD+OO 
-0.22232730-10 
0. OOOOOOOD+OO 


EPSO ( 1 ) - EPSO ( 6 ) 

O. 13989480-02 -0.41654050-04 

0 . OOOOOOOD +00 0 . OOOOOOOD + 00 


KAPP A ( 1 ) - KAPP A ( 6 ) 

-0 . 99 1042 ID- 1 7 0 . 5265 107D- 1 9 

0 . OOOOOOOD+OO 0 . OOOOOOOD+OO 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E(Y) E(Z) 


1 

2 

3 

4 


0 . 13989D-02 
0. 139890-02 
0. 13989D-02 
0. 13989D-02 


-0. 416540-04 
-0.416540-04 
-0.41654D-04 
-0.416540-04 


-0. 18252D-03 
-0. 182520-03 
-0. 18252D-03 
-0. 18252D-03 


-0.67079D- 1 1 
-0.67079D- 1 1 


S( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


5878379D-04 
. 67078 18D- 1 1 


. 1 294247D - 1 8 
. 4794883D-25 


E ( X Y ) 

-0.67078D- 1 1 
-0.67078D-1 1 
-0.67078D-11 
-0.67078D-1 1 
-0.67078D- 1 1 
-0 . 67078D- 1 1 


S ( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


. 18251810-03 
. 14819400-10 


. 57679510-18 
. 1 228763D-25 


E ( XY ) 

-0. 14819D-10 
-0. 14819D-10 
-0. 14819D-10 
-0. 1 48 1 90- 10 
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5 

0. 139890-02 

-0. 4 1654D-04 

-0. 18252D-03 

-0. 148 19D- 10 

6 

0. 139890-02 

-0.41654D-04 

-0. 18252D-03 

-0. 14819D-10 

THE STRESSES AT THE 

OUTER FIBER ARE 



PLY NO 

S( 1 ) 

5(2) 

S ( 3 ) 

S ( 6 ) 

1 

28287.59 

538 .68 

311.89 

0.00 

2 

-538.68 

1415.38 

311.89 

0.00 

3 

28287.59 

538.68 

311.89 

0.00 

4 

28287.59 

538.68 

311.89 

0.00 

5 

-538.68 

1415.38 

311.89 

0.00 

6 

28287 . 59 

538.68 

311.89 

0.00 

CYCLE 

NUMBER 0 

. 5999000E+04 




PLY NO. 


ALPHAM2 
28524400- 1 1 
55887680-03 
28524400- 1 1 
2852440D - 1 1 
5588768D -03 


0.28524400-1 1 


ALPHAM8 

0 . OOOOOOOD + 00 
-0. 2442928D- 10 
0 . OOOOOOOD +00 
0 . OOOOOOOD +00 
-0. 2442928D- 10 
0 . OOOOOOOD +00 


EPSO ( 1 ) - EPSO ( 6 ) 

0. 1402620D-02 
0 . OOOOOOOD +00 

KAPPA ( 1 ) - KAPP A ( 6 ) 


-0.4079366D-04 
0 . OOOOOOOD +00 


-0. 1831276D-03 
-0. 14295150-10 


-0 . 9590869D- 1 7 
0 . OOOOOOOD +00 


0. 5540902D- 19 
0 . OOOOOOOD + 00 


0 . 5245836D- 1 8 
0. 1 538825D -25 


THE STRAINS AT THE OUTER FIBER ARE 


NO . 

E ( X ) 

E ( Y ) 


E(Z) 

E(XY) 

1 

0. 14026D-02 

-0. 40794D-04 

-0. 

183 1 3D-03 

-0. 142950-10 

2 

0. 14026D-02 

-0.40794D-04 

-0. 

183130-03 

-O. 14295D- 10 

3 

0. 14026D-02 

-0. 407940-04 

-0. 

1 83 1 3D -03 

-0. 142950-10 

4 

0. 140260-02 

-0.40794D-04 

-0. 

1 83 1 3D -03 

-0. 14295D- 10 

5 

0. 14026D-02 

-0. 407940-04 

-o. 

183130-03 

-0. 142950-10 

6 

0 . 14026D-02 

-0.40794D-04 

-0. 

183 13D-03 

-0. 14295D- 10 

STRESSES AT THE 

OUTER FIBER ARE 




NO . 

S( 1 ) 

S ( 2 ) 


S( 3) 

S( 6 ) 

1 

28362.21 

541.66 

312 

.50 

0.00 

2 

-541 .66 

1340.76 

312 

.50 

0.00 

3 

28362.21 

541.66 

312 

.50 

0.00 

4 

28362.21 

541.66 

312 

.50 

0.00 

5 

-541 .66 

1340.76 

312 

.50 

0.00 

6 

28362 . 2 1 

541.66 

312 

.50 

0.00 


CYCLE NUMBER 


0 . 6499000E+04 


PLY NO. 


1 

2 

3 

4 

5 

6 


ALPHAM2 
0 . 4479322D- 1 1 
0.5870540D-03 
0.4479322D- 1 1 
0 . 4479322D- 1 1 
0.5870540D-03 
0 . 4479322D- 1 1 


ALPHAM8 

0 . OOOOOOOD +00 
-0. 2566094D-10 
0 . OOOOOOOD+OO 
0 . OOOOOOOD+OO 
-0 . 2566094D- 10 
0. OOOOOOOD+OO 


EP$0( 1 ) - EPSO ( 6 ) 


0. 1404679D-02 
0 . OOOOOOOD+OO 


-0.4031122D-04 
0 . OOOOOOOD+OO 


-0. 1834694D-03 
-0. 14001 19D-10 


KAPPA ( 1 ) - KAPPA ( 6 ) 


-0. 976 35 6 9D -17 
0 . 00000000+00 


0. 10967 10D- 18 
0. OOOOOOOD +00 


0 . 5457380D- 1 8 
0.4008539D-25 


THE STRAINS AT THE OUTER FIBER ARE 
PLY NO . E ( X ) E( Y) 


1 0.14047D-02 

2 0 . 1 4047D-02 

3 0. 140470-02 

4 0. 140470-02 


-0.4031 ID-04 
-0.4031 ID-04 
-0.4031 ID-04 
-0.4031 ID-04 


E(Z) 

-0. 18347D-03 
-0. 18347D-03 
-0. 183470-03 
-0. 18347D-03 


E(XY) 

-0. 1400 ID - 10 
-O. 1 400 1 D - 10 
-0. 14O01D-10 
-0. 1 4001 D- 10 
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5 O.14O47D-02 -0.4031 ID-04 -0.18347D-03 

6 O.14O47D-02 -0.4031 ID-04 -0. 183470-03 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO . S ( 1 ) S( 2 > S( 3 ) 


1 28404.05 

2 -543.33 

3 28404.05 

4 28404.05 

5 -543.33 

6 28404.05 


543 . 33 
1298 . 92 

543.33 

543 . 33 
1298.92 

543 . 33 


312.84 

-312.84 

312.84 

312.84 

-312.84 

312.84 


CYCLE NUMBER 


0 . 69990O0E+04 


PLY NO. 

1 

2 

3 

4 

5 

6 

EPSO ( 1 ) - E PSO ( 6 ) 

0. 1406102D-02 
O.OOOOOOOD+OO 


ALPHAM 8 

O.OOOOOOOD+OO 
-O. 26511670-10 
O.OOOOOOOD+OO 
0 . OOOOOOOD + 00 
-0. 2651 1670-10 
0 . OOOOOOOD +00 


-0. 3997799D-04 -0 

O.OOOOOOOD+OO -O 


ALPHAM2 
0.61 64068D - 1 1 
0. 60651630-03 
0.61 64068D - 1 1 
0 . 6 1 64068D - 1 1 
0.6065163D-03 
0 . 6 1 64068D - 1 1 


KAPPA! 1 ) - KAPP A ( 6 ) 

-0. 1000685D-16 0 . 1 477874D- 1 8 0 

O.OOOOOOOD+OO O.OOOOOOOD+OO 0 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E ( Y ) E(Z) 


1 

2 

3 

4 

5 

6 


O. 14061D-02 
0 . 1406 ID-02 
0. 1406 ID-02 
0. 1 406 1 D -02 
0. 1 406 1 D -02 
O. 1 406 1 D -02 


-0 39978D-04 
-0. 39978D-04 
-0. 39978D-04 
-0. 39978D-04 
-0.39978D-04 
-0- 399780-04 


-O. 183710-03 
-0. 1837 ID-03 
-0. 1837 ID-03 
-0. 1837 ID-03 
-0. 1837 ID -03 
-0. 18371D-03 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO . S ( 1 ) S(2) S( 3 ) 


1 28432 95 

2 -544.48 

3 28432.95 

4 28432.95 

5 -544.48 

6 28432.95 


544 . 48 
1270.02 

544 . 48 

544.48 
1270.02 
544 48 


313.08 

-313.08 

313.08 

313.08 

-313.08 

313.08 


CYCLE NUMBER 0 . 7499000E+04 


PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0 . 7892060D- 1 1 
0.62131 59D -03 
0. 7892060D-1 1 
O. 7892060D- 1 1 
0.62131 59D-03 
0 . 7892060D - 1 1 


ALPHAM 8 

O.OOOOOOOD+OO 
-0. 27 15858D- 10 
0 . OOOOOOOD +00 
0 . OOOOOOOD +00 
-0.2715858b- 10 
O.OOOOOOOD+OO 


EPSO ( 1 ) - EPS0(6) 



0. 1407183D-02 
0 . OOOOOOOD +00 

-0. 3972460D-04 -< 

0 . OOOOOOOD +00 -< 

KAPPA ( 1 ) - KAPPA ( 6 ) 




-0.9869425D-17 
0 . OOOOOOOD +00 

0 . 997492 ID- 19 < 

O.OOOOOOOD+OO i 

THE 

PLY 

STRAINS AT THE OUTER FIBER ARE 
NO . E ( X ) E ( Y ) 

E(Z) 


1 0 . 1 407 2D -02 

2 0 . 1 407 2D -02 

3 0. 14072D-02 

4 0.14072D-02 

-0. 39725D-04 
-0.39725D-04 
-0. 39725D-04 
-0.39725D-04 

-0. 18388D-03 
-0. 18388D-03 
-0. 18388D-03 
-0. 18388D-03 


-0. 14001D-10 
-0. 1 400 ID- 10 


$( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


. 1837055D-03 
. 1379814D-10 


. 5800803D-18 
. 567697 ID-25 


E(XY) 

-0. 13798D- 10 
-0. 13798D- 10 
-0. 13798D-10 
-0. 13798D-10 
-0. 13798D- 10 
-0. 13798D- 10 


$( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


. 1838850D-03 
. 1 364374D- 10 


. 5631240D- 18 
. 3486876D-25 


E(XY) 

-0. 13644D-10 
-0. 13644D-10 
-0. 1 3644D- 10 
-0. 13644D-10 
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5 

0. 1 407 2D -02 

-0.39725D-04 -0- 18388D-03 

-0. 13644D-10 


6 

0. 14072D-02 

-0.39725D-04 -0. 18388D-03 

-0. 136440-10 

THE 

STRESSES AT THE 1 

3UTER FIBER ARE 


PLY 

NO. 

S( 1 ) 

512) S ( 3 ) 

S( 6 ) 


1 

28454.92 

545.35 313.26 

0.00 


2 

- 545 . 35 

1248.05 -313.26 

0.00 


3 

28454.92 

545.35 313.26 

0.00 


4 

28454.92 

545.35 313.26 

0.00 


5 

-545.35 

1248.05 -313.26 

0.00 


6 

28454.92 

545.35 313.26 

0.00 

CYCLE NUMBER 0 

. 7999000E +04 


PLY 

NO . 

ALPHAM2 ALPHAM8 



1 

0.985463 ID- 1 1 0 . OOOOOOOD + 00 



2 

0.6332 189D -03 -0. 2767888D- 10 



3 

0.965463 ID- 1 1 0 . OOOOOOOD+OO 



4 

0.965463 ID- 1 1 0 . OOOOOOOD+OO 



5 

0.6332 189D -03 -0 . 2767888D - 1 0 



6 

0. 965463 ID- 1 1 0 . OOOOOOOD+OO 


EPSO( 1 ) 

- EPSO ( 6 ) 




0 

. 1 408053D -02 

-0.39520800-04 -0 

. 1840293D-03 


0 

. OOOOOOOD+OO 

0. OOOOOOOD+OO -0 

. 1351956D-10 

KAPPA ( 1 ) 

- K APPA (6 ) 




-0 

, 9686882D- 17 

0 . 6966505D - 1 9 0 

.53705260-18 


0 

. OOOOOOOD+OO 

0 OOOOOOOD+OO 0 

! 2 1681 40D-25 

THE 

STRAINS AT THE OUTER FIBER ARE 


PLY 

NO . 

E ( X ) 

E(Y) E(Z) 

E(XY ) 


1 

0. 1 408 ID -02 

-0.39521D-04 -0. 18403D-03 

-0. 13520D-10 


2 

0. 1 408 1D~02 

-0 3952 1 D -04 -0. 18403D-03 

-0. 13520D- 10 


3 

0. 1 408 ID-02 

-0. 395210-04 -0. 18403D-03 

-0. 13520D- 10 


4 

0. 1 408 ID -02 

-0. 3952 ID -04 -0. 18403D-03 

-0. 13520D-10 


5 

0, 1 408 ID-02 

-0 . 3952 ID-04 -0. 18403D-03 

-0. 13520D-10 


6 

0. 1408 ID -02 

-0 . 3952 ID-04 -0. 18403D-03 

-0. 1 3520D- 10 

THE 

STRESSES AT THE OUTER FIBER ARE 


PLY 

NO . 

S( 1 ) 

S ( 2 ) S ( 3 ) 

S(6 ) 


1 

28472.60 

546.06 313.40 

0.00 


2 

-546.06 

1230,37 -313.40 

0.00 


3 

28472.60 

546.06 313.40 

0.00 


4 

28472.60 

546.06 313.40 

0.00 


5 

-546.06 

1230.37 -313.40 

0.00 


6 

28472.60 

546.06 313.40 

0.00 

CYCLE NUMBER 0 

. 8499000E+04 


PLY 

NO . 

ALPHAM2 ALPHAM8 



1 

0. 1 1 44603D - 1 0 0. OOOOOOOD+OO 



2 

0. 64315 17D -03 -0.281 1306D-10 



3 

0. 1 144603D-10 0. OOOOOOOD+OO 



4 

0. 1 144603D- 10 0. OOOOOOOD+OO 



5 

0. 6431517D-03 -0.281 1306D- 10 



6 

0. 1 1 44603D- 10 0. OOOOOOOD+OO 


EPSO( 1 ) 

- EPS0(6 ) 




0 

. 14087790-02 

-0. 39350730-04 -0. 

184 1498D-03 


0 

. 00000000+00 

0. OOOOOOOD+OO -o. 

1341594D-10 

KAPPA ( 1 ) 

- KAPPA(6 ) 




-o 

i 9947443D-17 

0.821 8032D- 1 9 0. 

5771 187D-18 


0 

' . OOOOOOOD+OO 

0. OOOOOOOD+OO 0. 

26 1 8052D-25 

THE 

STRAINS AT THE OUTER FIBER ARE 


PLY 

NO . 

E ( X ) 

E ( Y ) E ( Z ) 

E ( XY ) 


1 

0. 1408QD-02 

-0 . 3935 ID-04 -0.18415D-03 

-0. 13416D- 10 


2 

0. 140880-02 

-0.3935 ID-04 " -0. 18415D-03 

-0. 13416D-10 


3 

0. 14088D-02 

-0. 3935 ID-04 -0.18415D-03 

-0. 13416D-10 


4 

0. 14088D-02 

-0. 3935 ID-04 -0.18415D-03 

-0. 13416D-10 
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5 0 . 1 4088D-02 -0.393510-04 -0-184150-03 -0.13416D-10 

6 0.14088D-02 -0.39351D-04 -0.18415D-03 -0. 134160-10 

THE STRESSES AT THE OUTER FIBER ARE 

PLY NO. S ( 1 ) S ( 2 ) S ( 3 ) S(6) 

1 28487.35 546.65 313.52 0.00 

2 -546.65 1215.62 -313.52 0.00 

3 28487.35 546.65 313.52 0.00 

4 28487.35 546.65 313.52 0.00 

5 -546.65 1215.62 -313.52 0.00 

6 28487.35 546.65 313.52 0.00 

CYCLE NUMBER 0 . 8999000E+04 

PLY NO. ALPHAM2 ALPHAM8 

1 0 . 13262 14D- 10 O . 00000000+00 

2 0.65 1660 1 D -03 -0.28484970-10 

3 0.13262140-10 0 OOOOOOOD + 00 

4 0. 13262 14D- 10 0 . OOOOOOOD + 00 

5 0. 65166010-03 -0 . 2848497D- 10 

6 O. 13262 14D- 10 0 . OOOOOOOD+OO 

E PSO ( 1 ) - EPSO ( 6 ) 

0. 14094010-02 -0. 3920505D-04 -O. 1842530D-03 

0 . OOOOOOOD +00 0.00000000+00 -o. 1 3327 170-10 

KAPPA( 1 ) - KAPPA ( 6 ) 

-0 . 982434 ID- 1 7 0.70201570-19 0.55862030-18 

O . 00000000+00 0 . OOOOOOOD + 00 0.21 20855D-25 

THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E(X) E ( Y ) EfZ) E(XY) 

1 0.14094D-02 -0. 392050-04 -0.184250-03 -0. 133270-10 

2 0.14094D-02 -0. 392050-04 -0.184250-03 -0.13327D-10 

3 0.14094D-02 -0. 392050-04 -0.18425D-03 -0.133270-10 

4 0.14094D-02 -Q.39205D-04 -0.18425D-03 -0.13327D-10 

5 0.140940-02 -Q.39205D-04 -Q.18425D-03 -0.13327D-10 

6 0.140940-02 -0.39205D-04 -0.184250-03 -0.133270-10 

THE STRESSES AT THE OUTER FIBER ARE 
PLY NO . S( 1 ) S(2) S( 3) 

1 28499.98 . 547.15 313.62 

2 -547.15 1202.99 -313.62 

3 28499.98 547.15 313.62 

4 28499.98 547.15 313.62 

5 -547.15 1202.99 -313.62 

6 28499.98 547. 15 313.62 

CYCLE NUMBER 0 . 9499000E+04 

ALPHAM8 

0 . OOOOOOOD +00 
-0. 2880982D- 10 
0 . OOOOOOOD +00 
0 . 00000000+00 
-0.28809820-10 
0 . OOOOOOOD +00 

EPSO ( 1 ) - EPSO ( 6 ) 

0. 14099440-02 -0 . 390778 1 D-Q4 -0. 18434320-03 

0 OOOOOOOD +00 0 . OOOOOOOD +00 -O. 1324964D- 10 

KAPPA( 1 ) - KAPPA ( 6 ) 

-0 98336 15D- 17 0 . 78861930- 19 0.55902460-18 

0.00000000+00 0.00000000+00 0.25253330-25 

THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E(Y) E ( Z ) E(XY) 

1 0.140990-02 -0.390780-04 -0.184340-03 -0.13250D-10 

2 0.14099D-02 -O. 390780-04 -0.18434D-03 -013250D-10 

3 0.14099D-02 -0.39078D-04 -0.18434D-03 -0.13250D-10 

4 0.140990-02 -0.39078D-04 -Q.18434D-03 -0.13250D-10 


PLY NO. ALPHAM2 

1 0. 1509990D- 10 

2 0 . 65909 1 7D-03 

3 0.15099900-10 

4 0.15099900-10 

5 0.65909170-03 

6 0. 1509990D- 10 


S ( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 
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5 0.14099D-02 -O. 390780-04 -0.18434D-03 

6 0 140990-02 -0.39078D-04 -0.18434D-03 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO. 5(1) 5(2) SO) 


1 28511.02 

2 -547.59 

3 28511.02 

4 28511.02 

5 -547.59 

6 28511.02 


547.59 
1191.95 
547.59 
547.59 
1 191 .95 
547.59 


313.71 

- 313.71 

313.71 

313.71 
- 313.71 

313.71 


CYCLE NUMBER 


0 . 9999000E+04 


PLY NO. 

1 

2 

3 

4 

5 

6 


ALPHAM2 
0. 1695690D- 10 
0 . 6656820D-03 
0. 1695690D- 10 
0. 1 695690D- 10 
0. 66568200-03 
0. 1 695690D- 10 


ALPHAM8 — - - 

0 . 0000000 D +00 
-0. 2909789D- 10 
0 . OOOOOOOD + 00 
0 . OOOOOOOD + 00 
-o. 2909789D- 10 
0 . OOOOOOOD + 00 


EPSO ( 1 ) - EPSO( 6 ) 


0. 1410425D-02 
0 . OOOOOOOD +00 

KAPPA ( 1 ) - KAPPA ( 6 ) 

-0. 1 006 5 2 1 D - 1 6 
0 . OOOOOOOD +00 


-0.3896497D-04 -0 

0 . 00000000+00 -0 

-0. 1896032D-19 0 

O.OOOOOOOD+OO -0 


THE STRAINS AT THE OUTER FIBER ARE 

PLY NO. E ( X ) E ( Y ) E(Z) 


1 

2 

3 

4 

5 

6 


0. 141040-02 
0. 14104D-02 
0. 14 104D-02 
0. 14 104D-02 
0. 14 104D-02 
0. 14 104D-02 


-0.38965D-04 
-0. 38965D-04 
-0.38965D-04 
-0, 38965D-04 
-0.38965D-04 
-0. 38965D-04 


-0. 18442D-03 
-0. 18442D-03 
-0. 18442D-03 
-0. 18442D-03 
-0. 18442D-03 
-0. 18442D-03 


THE STRESSES AT THE OUTER FIBER ARE 

PLY NO. S( 1 ) S ( 2 ) SO) 


1 28520.81 

2 -547.98 

3 28520 81 

4 28520.81 

5 -547.98 

6 28520.81 


547.98 
1182.17 
547.98 
547.98 
1 182 . 17 
547 .98 


313.79 
-313.79 
3 13.79 
313.79 
-313.79 
313.79 


-0. 13250D- 10 
-0 1 3250D- 10 


S ( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 


. 184423 ID-03 
. 1318089D- 10 


. 6065785D- 1 8 
. 221561 2D- 25 


E(XY) 

-0. 1318 ID- 10 
-0. 13181D-10 
-0. 13181D-10 
-0. 1318 ID- 10 
-0. 1318 ID- 10 
-0. 1 3 18 ID- 10 


S( 6 ) 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 
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